ChinaXiv 合 作 期 刊 


Eb 和 a“ 双 追踪 ”GRB131231A 的 能 谱 特 性 研究 ' 


王 展 ， 彭 朝阳 +*， 陈 家 明 
(云南 师范 大 学 物理 与 电子 信息 学 院 ， 云 南 昆明 650500) 
摘要 : GRB 131231A 是 一 个 具有 友和 a“ 双 追踪 ”行为 的 特殊 暴 。 为 了 进一步 分 析 该 
暴 的 能 谱 特 征 ， 我 们 使 用 贝 叶 斯 方法 对 它 做 了 详细 的 时 间 分 辨 光谱 分 析 后 发 现 24 个 分 辨 谱 
中 16 个 有 明显 热 成 分 ， 这 些 热 成 分 主要 分 布 在 脉冲 的 前 期 和 峰值 附近 ， 支 持 了 伽 玛 暴 喷 流 
从 物质 主导 到 磁 主 导 的 过 渡 。 加 入 热 成 分 后 参数 Ep 、a 都 变 硬 ， 而 峰值 流量 F, 几乎 不 变 。 我 
们 使 用 光 球 层 成 分 计算 了 相对 论 流出 物 的 特性 参数 后 发 现 流出 物 特性 参数 初始 半径 mm 、 饱 和 
FREE rg. HRE rn 流出 物 滑 行 运动 阶段 洛 伦 兹 因子 r 均 发 生 较为 明显 的 演化 ; 其 温 
度 随时 间 的 演化 也 能 用 拐 折 窜 律 模型 拟 合 ， 其 上 升 和 下 降 阶段 的 究 律 指数 分 别 为 a = 7. 65 
土 3. 48、b = -2. 10 土 3. 48， 下 降 指 数 b 与 前 人 得 出 的 接近 于 -2/3 有 较 大 偏差 。 
关键 词 : Fermi 伽 玛 射线 暴 ， 瞬 时 辐射 能 谱 ， 光 谱 拟 合 ， 数 据 分 析 ， 最 优 模 型 
中 图 分 类 号 : P172.3 文献 标识 码 : A 文章 编号 : 
0. 引 言 
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. 依然 让 人 感到 困惑 , 一 些 非常 基础 的 问题 是 喷 流 成 分 (重子 主导 的 火球 或 是 坡 印 廷 流量 主导 
ia 的 流出 物 )， 它 决定 或 至 少 强烈 影响 能 量 耗 散 机 制 〈 激 波 或 磁 重 联 )， 粒 子 加 速 机 制 〈 热 主导 
或 磁 主 导 )， 特 别 是 辐射 机 制 〈 同 步 辐射 ) 或 光 球 准 热 辐射 的 康 普 顿 化 品 。 伽 玛 暴 辐射 有 瞬 
时 辐射 和 余辉 两 个 主要 阶段 。 瞬 时 辐射 阶段 通常 持续 几 秒 到 几 分 钟 。 伽 玛 暴 瞬时 辐射 具有 和 丰 
富 的 观测 特征 , 但 是 其 中 的 大 多 数 没 有 被 很 好 的 解释 。 伽 玛 暴 的 瞬时 辐射 阶段 的 光谱 成 分 主 
要 来 自 两 种 起 源 : 同步 辐射 成 分 〈 非 热 成 分 ) 和光 球 层 ( 热 成 分 ) 叫 。 过 去 的 研究 表明 ， 
GRBs 能 谱 中 可 能 存在 三 种 基础 成 分 申 : (1) 能 被 Band 函数 拟 合 的 非 热 成 分 (2) 能 被 普 
朗 克 函数 拟 合 的 热 成 分 ;3) 出 现在 能 量 较 高 区 域 的 震 律 成 分 。 目 前 已 观测 到 的 GRBs f 
够 用 发 现 两 种 或 两 种 以 上 的 成 分 外 。 热 成 分 主导 的 GRBs 非常 稀少 ， 所 以 对 于 暴 的 能 谱 处 理 
时 多 采用 组 合 模型 ， 常 用 的 基础 模型 通常 是 Band 或 者 CPL。 多 数 情 况 下 , 组 合 模型 往往 
得 到 更 优 的 参数 值 ， 这 更 有 助 于 我 们 推断 GRB 能 源 机 制 和 中 心 引擎 。 在 BATSE 时 期 ，Ep 
显示 出 几 种 不 同 的 演化 : (1)“ 从 硬 到 软 ” 的 趋势 ， 无 论 Flux 的 升降 如 何 ， 它 都 是 单调 递减 
的 PMS, (2) “流量 跟踪 ”趋势 ， 即 参数 的 演化 跟随 Flux 上 升 和 下 降 中 。(3)“ 软 到 硬 ” 的 趋势 
或 混乱 的 演变 ) 钙 。2008 年 费 米 卫星 升 空 后 ， 大 量 的 伽 玛 暴 观测 数据 证 实 前 两 种 模式 占 主导 
地 位 :由 硬 到 软 ” 约 占 三 分 之 二 ，“Flux 跟踪 ” 约 占 三 分 之 一 站。 这 些 演化 模式 的 物理 起 源 仍 
然 没 有 解决 ， 尽 管 文献 中 已 经 提出 了 一 些 方案 MO, EG Ep 演化 相 比 ，w 演化 更 加 混乱 ， 
因此 研究 和 物理 解释 相对 较 少 。 李 亮 发 现 GRB 131231A 的 Ep 和 演化 都 表现 出 “ 通 量 跟 
踪 ” 行 为 ， 并 且 定 义 为 光谱 演化 的 “ 双 跟 踪 ” 模 式 中 。 但 李 亮 等 并 没有 加 入 热 成 分 进行 详 
的 能 谱 分 析 , 既然 很 多 研究 表明 伽 玛 暴 能 谱 中 包含 热 成 分 , 为 了 更 详细 地 研究 光谱 参数 Ep 和 
a 的 “ 双 追 踪 ” 情 况 的 能 谱 特 性 ， 我 们 使 用 贝 叶 斯 分 区 并 使 用 MCMC (蒙特 卡 罗 模 拟 方 法 ) 
处 理 GRB131231A。 在 GRB 的 研究 中 , 研究 人 员 发 现 次 主导 的 热 成 分 ， 一般 出 现在 Band 谱 
0 的 左 侧 03， 即 热 成 分 一 般 出 现在 低能 Band 谱 区 域 。 针 对 这 种 现象 ， 李 亮 在 探究 八 个 具有 
较 明 显 热 成 分 的 GRBs 中 热 成 分 对 于 非 热 成 分 中 几 个 特征 参数 值 的 影响 时 ， 都 选用 CPL 作 
为 基础 模型 与 BB 模型 组 合 来 进行 探究 巴 。 对 于 绝 大 多 数 伽 玛 暴 的 时 间 分 辨 谱 ，CPL 模型 是 
适合 的 5505。 然 而 在 一 些 例 子 中 ，Band 模型 的 高 能 寡 律 指数 B 显著 地 提高 拟 合 。 在 同步 辐 
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射 的 情况 下 ， 高 能 露 律 指数 B 提供 了 有 关 粒 子 加 速 和 激 波 的 本 质 ， 而 在 光 球 辐射 的 情况 下 ， 
它 提供 了 有 关 辐 射 区 域 能 量 耗 散 的 信息 59。 考 虑 到 这 种 状况 ， 我 们 以 统计 值 BIC 作为 判断 
标准 5， 对 伽 玛 暴 的 每 个 拟 合 区 间 选 用 最 优 模 型 的 数据 ， 并 研究 使 用 最 优 模 型 之 后 
GRB131231A 中 能 谱 参 数 的 变化 。 


1 数据 处 理 和 分 析 方 法 
GRB131231A( 触 发 号 : 131231198) 由 Fermi 卫星 中 的 伽 玛 暴 探测 器 (GBM: 8 keV - 40 MeV) 
于 2013 的 12 月 31 触 发 的 。 本 工作 中 ,我 们 使 用 贝 叶 斯 分 析 包 ， 即 最 大 似 然 工作 框架 (3ML) 
来 对 所 有 伽 玛 暴 的 光谱 分 析 。 该 工具 包 是 采用 马尔 可 夫 链 蒙特 卡 罗 理 论 (Markov Chain Monte 
Caro，MCMC) 对 时 间 分 辩 谱 进行 处 理 , 该 种 方法 的 优点 请 参阅 鼎 。 我 们 获取 费 米 卫星 的 
GBML3 来 自 网 站 (https:/fermi.gsfc.nasa.gov/)， 在 GBM 的 数据 中 有 三 种 数据 类 型 CTIME 
文件 、CSPEC 文件 和 TTE 文件 ， 我 们 选择 具有 最 小 时 间 分 辨 率 (2ns) 的 TTE 文件 来 进行 妆 
据 分 析 。 在 该 贝 叶 斯 统计 的 MCMC 拟 合 例子 中 ， 为 了 确定 能 谱 拟 合 的 最 优 模型 ， 我 们 选取 
BIC (Bayesian information criterion)! 作为 判断 标准 ， 我 们 通过 比较 不 同 模型 BIC 的 差 值 
CABIC) 来 确定 最 佳 模型 。 一 般 来 说 ，Band 模型 ”和 CPL 模型 ”在 伽 玛 暴 拟 合 时 是 两 种 
较为 优越 的 模型 ， 我 们 把 这 两 种 模型 选 作 基础 模型 ， 让 它们 分 别 与 BB RBA, RATED 
玛 暴 的 时 间 分 辨 谱 的 每 一 个 区 间 都 经 过 模型 对 比 , 进而 选择 出 该 区 间 适 用 的 最 优 模 型 。 为 了 
更 好 地 研究 光谱 随时 间 的 演化 , 对 光谱 进行 时 间 分 辨 谱 分 区 是 一 种 有 效 的 研究 方法 ,Burgess 
研究 过 四 种 对 光 变 曲 线 分 时 间 片 的 工具 9， 其 中 信 噪 比 (S/N) 和 贝 叶 斯 分 区 (BBs) "是 两 种 
最 有 效 的 分 区 方法 , 他 指出 这 两 种 方法 各 有 优 缺 点 ,传统 的 信 噪 比分 区 可 以 确保 每 个 时 间 段 
能 分 到 足够 多 的 光子 去 做 光谱 拟 合 , 但 是 有 时 它 会 破坏 物理 结构 。 关 于 贝 叶 斯 分 区 的 优 缺 点 
请 参阅 [19]， 我 们 采用 贝 叶 斯 区 间 GBBs) 分 区 ， 但 是 依然 使 用 信 噪 比 作为 一 种 筛选 方法 去 易 
除 所 得 区 间 中 不 满足 某 一 标准 的 区 间 (本 工作 中 我 们 选取 信 噪 比 S/N 三 20 的 区 间 ， 因 为 信 噪 
比 S/N<20 的 区 间 在 拟 合 时 会 出 现 较 大 的 误差 品 , 这 两 种 方法 的 共同 使 用 让 我 们 得 到 更 有 研 
究 意 义 的 区 间 。 表 1 给 出 了 该 暴 的 相关 参数 , 注 : CPL! 代表 ABIC<-8 的 时 间 分 辨 谱 区 间 使 
用 CPL #22, CPL? 代表 -8<ABIC<8 的 时 间 分 辨 谱 区 间 使 用 CPL 模型 。 图 1 给 出 了 一 个 例 
子 分 辨 谱 拟 合 结果 。 


表 1 件 玛 暴 参 数 及 使 用 的 探测 器 编号 
Table.1 Gamma-ray burst parameters and detector numbers used 


有 热 成 分 EA CPL' 区 间 EA CPL’ 区 间 数 ”适用 Band 区 间 
的 区 间 数 。” BC ABIC<-8) (-8< ABIC<8) 数 (8<ABIC) 
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1 一 个 黑体 成 分 最 显著 的 计数 率 时 间 分 辨 谱 (时间 区 间 为 24.59-25.86) 拟 合 结 
R ERN Band 模型 ， 右 图 为 Band+bb) 
Fig.1 Example of count spectral fit results with the most significant black body 
components in the time bin 24.59s-25.86s, where the left panel is band model, and the 


chinaXiv:202204.00113v1 


3. 分 析 结 果 


我 们 使 用 的 波段 函数 模型 、 有 局 
拟 合 ， 得 到 了 它们 的 时 间 分 罗 
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right panel is Band+bb. 


Rin ee A. PEA oe ee A, X} GRB131231A 进行 
eR, UZ 2。 我 们 从 时 间 分 辨 谱 的 ABIC 来 看 ， 有 些 暴 的 


某 些 时 间 段 中 CPL 模型 的 拟 合 结果 优 于 Band 模型 ， 有 些 结果 相反 ; 综合 考虑 之 后 我 们 把 每 


个 暴 的 每 个 时 间 分 辩 谱 区 间 进 行 最 优 模型 选择 ， 我 们 所 用 的 选择 方法 是 用 每 个 区 间 单 独 用 


CPL 模型 拟 合 时 产生 的 BIC 减 去 单独 用 Band 模型 拟 合 时 产生 的 ABIC 作为 筛选 标准 ， 当 
ABICS-8 时 我 们 认为 该 区 间 更 适合 CPL 模型 ， 当 ABIC>8 时 我 们 认为 该 区 间 更 适合 Band 
BAA, ~4-8<ABIC<8 时 我 们 认为 该 区 间 更 适合 CPL 模型 ( 当 两 种 模型 中 任何 一 种 模型 优势 不 


明显 时 ， 一 般 是 选用 参数 少 的 模型 ) 。 然 后 再 使 用 该 区 间 选 定 基础 模型 的 BIC 与 基础 模型 加 


上 BB 模型 的 组 合 模型 的 BIC 差 值 是 否 大 于 等 于 8 来 判断 该 区 间 是 否 有 热 成 分 , 最 后 筛选 出 


信 品 比 SIN S20 的 区 


光谱 分 析 时 , 我 们 
每 个 时 间 区 间 中 参数 的 变化 是 准 亚 
区 间 参 数 的 行为 号 。 在 文中 ， 我 们 展示 了 伽 玛 暴 的 


3.1 参数 的 分 布 


间 ， 作 为 研究 对 象 , WK 4。 


得 到 的 单个 时 间 区 间 的 时 间 跨 度 远 小 于 整体 的 时 间 跨 度 , 我 们 假设 在 


态 的 , 我 们 用 时 间 区 间 中 值 时 刻 参 量 的 特征 代表 整个 时 间 


能 谱 参数 演化 图 。 


图 2 显示 了 Best 模型 参数 分 布 ， 包 括 aq、Ep、Fp。 表 2 列 出 了 分 布 的 平均 值 和 标准 偏 
差 。 我 们 对 各 个 参数 分 布 进行 了 核 密度 估计 (KDE), 


入 黑体 成 分 和 加 入 了 黑体 成 分 时 低 


采用 高 斯 函数 核 。 如 表 2 所 示 ， 未 加 


能 光子 指数 w 分 布 的 均值 分 别 为 :-1.06+0.29,-1.00+0.29， 


加 入 黑体 成 分 后 ，a 轻微 地 变 硬 了 ， 但 他 们 都 与 典型 值 ~ -1 基本 一 致 ， 且 都 没有 超过 所 
谓 的 同步 “ 死 线 ”(-3/2<a<-2/3)。 在 加 入 黑体 后 ， 峰 值 能 量 5p 的 核 密度 曲线 的 峰值 向 右 发 生 
了 偏 移 ， 均 值 为 199.84 土 161.61， 和 未 加 入 黑体 成 分 相 比 ， 峰 值 能 量 E, 确实 也 变 硬 了 ; 而 
F, 加 入 黑体 前 后 没有 明显 变化 ， 平 均值 仅 有 微弱 变化 ， 这 些 结果 和 前 人 的 研究 结果 一 致 。 
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(b) 
图 2 Best 模型 的 参数 a、Ep、Fp 分 布 
Fig.2 The distribution of the parameter values Ep , a and Fp. 


表 2 模型 参数 均值 与 标准 偏差 


Table 2 Mean and standard deviation of model parameters. 
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3.2 MARA AIS 


普 参 数 的 演化 
我 们 用 我 们 分 析 的 数据 重新 分 析 未 加 入 黑 


= 
— 
= 


的 峰值 能 量 和 低能 谱 指数 的 演化 情况 ， 
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图 3 中 展示 贝 叶 斯 区 间 中 S$ 三 20 的 时 间 分 辩 谱 区 间 对 应 的 数据 点 .GRB131231A 的 参数 Eb 和 
co 表 现 为 较 好 的 “流量 跟踪 ”趋势 ， 我 们 得 到 的 结果 和 李 亮 的 结果 一 致 。 
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图 3. GRB131231A 的 最 优 模型 的 时 间 分 辨 谱 的 参数 值 5p 值 和 a 值 的 分 布 。 
Fig.3 The distribution of the parameter values Ep and a of the time-resolved spectrum of 
the optimal model of GRB131231A. 


3. 3 加 入 热 成 分 后 光谱 参数 的 演化 及 变化 

当 加 入 黑体 成 分 后 ， 我 们 找 出 探测 到 热 成 分 的 16 个 时 间 分 辩 谱 区 间 发 现 ， 两 个 参数 Ep 
Al 的 演化 趋势 和 未 加 入 热 成 分 的 趋势 一 致 ， 依 然 表 现 出 双 追 踪 行 为 。 如 3. 1 节 所 述 值 ， 这 
两 个 参数 发 生 了 变化 ，a 和 EE 轻微 变 硬 ， 但 从 图 4 中 我 们 能 发 现 这 两 个 参数 变化 的 细节 ， 
GRB131231A 中 有 热 成 分 的 区 间 的 Ep 值 的 整体 变化 并 不 一 致 ， 如 图 4(a) 。 我 们 发 现 
GRB131231A 中 Ep 值 在 峰值 之 前 ， 加 入 BB 成 分 后 变 硬 〈 值 变 大 ); 峰值 之 后 ， 加 入 BB 成 
分 后 变 软 〈 值 变 小 )。 其 中 Ep 值 变 硬 的 点 有 6 个 ， 均 分 布 在 峰值 之 前 ，Ep 值 变 软 的 点 有 10 
个 ， 均 分 布 在 峰值 后 。GRB131231A 中 wa 的 整体 变化 也 不 一 致 ， 如 图 4 (b)。 有 一 个 点 分 布 
在 同步 辐射 死 线 之 上 其 余 的 点 都 分 布 在 同步 辐射 死 线 之 下 ; 在 加 入 BB 成 分 后 ，a 出 现 了 明 
显 的 变化 ， 在 峰值 之 前 ， 加 入 BB 成 分 后 a 值 变 软 ;， 峰值 之 后 ， 加 入 BB 成 分 后 变 硬 。 在 加 
入 黑体 前 ，GRB131231A 中 仅 有 最 后 两 个 点 的 ax 值 小 于 - 3/2; 加 入 热 成 分 之 后 ， 脉 冲 峰值 
处 有 两 个 点 的 a 值 上 升 到 超过 同步 辐射 死 线 ， 最 后 两 个 点 的 wx 值 略 有 上 升 ， 但 是 依然 小 于 
一 3/2。 
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图 4. GRB131231A 中 探测 到 热 成 分 的 时 间 分 辨 谱 区 间 中 Ep 值 和 a 值 的 分 布 ,(b) 图 上 下 两 
条 红色 虚线 分 别 代表 a= 一 2/3 和 a= 一 3/2 的 位 置 。 

Fig.4 Distribution of Ep and a in the time-resolved spectral interval of thermal components 
detected in GRB131231A, the upper and lower red dashed lines represent o=-2/3 and a=- 3/2. 


表 3 GRB131231A 在 不 同 模型 中 Ep-F,。、Ep-q、Fp-a 人 参数 关系 
(Spearman 相关 性 系数 r+ 和 置信 和 度 p) 
Table 3 GRB131231A parameter relationship of Ep-Fy, Ep-a and 
F,-a in different models 


(Spearman correlation coefficient r and confidence p) 


模型 Ep-Fp Ep-a Fp-a 
p r p r p 

Best 0.426 0.0378 0.412 0.0453 0.765 < 107+ 
Best +BB 0.332 0.208 0.065 0.812 0.738 0.001 

Band 0.687 0.00021 0.70 0.00014 0.655 0.00052 
Band +BB 0.444 0.085 0.476 0.062 0.871 < 1074 

CPL 0.566 0.0039 0.731 < 1074 0.756 < 10+ 
CPL+BB 0.368 0.161 -0.082 0.762 0.70 0.003 


3. 4 光谱 参数 的 相关 性 

相关 分 析 在 理解 GRB 物理 学 中 起 着 重要 作用 ， 因 为 它 提 供 了 一 些 线索 来 帮助 我 们 揭示 
它们 的 本 质 。 我 们 即 log(E,)-log(F,)、log(F,)-a、log(E,)-a， 表 3 具体 给 出 了 相关 系数 。 相 
关 的 强度 由 Spearman 的 系数 r 给 出 ， 强 相关 性 r> 0.7， 弱 相关 性 r <0.4。 [7] [24] [25] 揭 示 了 
log(F)-a 关系 主要 有 三 种 类 型 : 非 单调 关系 (包含 正和 负 寡 律 段 ， 在 峰值 流量 处 典型 地 具有 了 明 
显 的 间断 )、 单 调 关 系 (由 单个 暴 律 描述 ) 和 没有 明显 的 趋势 。GRB131231A 三 种 模型 中 
log(F,)-log(E,) 相 关 性 Best: r=0.426, p=0.0378 , Band:r=0.687, p<10*, CPL: r=0.566, p<10*, 
最 好 的 是 Band 模型 ， 三 种 模型 中 log(E,)-a 相关 性 Best: r=0.412, p=0.0453, Band:r=0.70， 
p=0.00014, CPL: r=0.731，p<10“ 最 好 的 是 Band 模型 ， 三 种 模型 中 log(F)-a 相关 性 Best: 
r=0.765, p<10*, Band:r=0.655, p=0.00052, CPL: r=0.756，p<10“ 最 好 的 是 Best 模型 。 在 
Band 和 CPL 模型 中 log(E,)-log(F), log(F,)-a. log(E,)—o 都 表现 出 极 强 的 相关 性 ， 当 我 们 使 
] Best 模型 时 ,log(F,)-log(E,)、log(E,)-o 的 相关 性 均 变 差 ;Best 模型 的 log(F,)-a 相关 性 强 于 
Band 和 CPL。 
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图 5 (a) ABIC 随 脉冲 的 演化 和 分 布 图 , 绿色 和 红色 虚线 分 别 表 示 ABIC 等 于 0、-8 的 位 置 ; (b) 
探测 到 热 成 分 的 光谱 中 热 成 分 流量 与 总 能 量 流量 的 比率 ， 红 色 虚 线 代 表 比 率 为 1。 
Fig.5 (a) The evolution and distribution of A BIC with pulse, and the green and red dotted lines 


represent the positions where A BIC is equal to 0 and - 8 respectively; (b) The ratio of thermal 
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component flow to total energy flow in the spectrum of detected thermal component, and the red 


dotted line represents the ratio of 1. 


3.5 热 成 分 分 析 

为 了 检查 能 谱 中 的 热 成 分 ， 我 们 先 寻 找 对 应 于 时 间 分 辨 谱 的 最 优 模 型 ， 接 着 用 最 优 模 
型 加 入 BB 成 分 拟 合 ， 查 看 加 入 黑体 前 后 统计 量 (BIC) 的 差 值 ABIC， 如 果 ABIC <-8 说 明光 
谱 中 显著 含有 热 成 分 。 图 5 中 a 图 给 出 了 ABIC 在 脉冲 中 的 演化 和 分 布 图 ，b 图 给 出 了 热 成 
分 流量 与 总 能 量 流量 的 比率 。 

GRB131231A 具有 明显 的 脉冲 上 升 阶 段 和 脉冲 下 降 阶 段 ， 我 们 对 GRB131231A 中 出 现 
热 成 分 的 区 域 进行 了 统计 ， 从 a 图 统计 的 结果 来 看 ABIC <-8 的 光谱 主要 出 现在 对 应 高 的 光 
子 计 数 的 光 变 曲线 的 位 置 ， 这 说 明 热 成 分 主要 存在 于 脉冲 前 期 ， 并 且 随 着 时 间 逐 渐 减 少 。b 
图 表明 热 成 分 流量 与 总 能 量 流量 的 比率 多 数 分 布 在 0. 1 之 下 , 比率 较 高 的 区 域 主要 分 布 在 高 
脉冲 区 域 。 


3. 6 与 热 成 分 相关 的 参数 值 
> 现在 的 研究 表明 热 成 分 起 源 于 光 球 层 , 根据 高 纬 辐射 解释 , 最 高 温度 和 最 强 热流 最 初 从 
CD 光 球 层 辐射 轴 上 发 出 ， 并 朝向 观测 者 。 因 此 从 原则 上 来 说 ， 光 子 辐射 位 置 的 半径 可 以 被 确定 
P9。 按 照 [26] 的 理论 ， 对 于 一 个 已 知 红 移 的 伽 玛 暴 当 我 们 知道 了 该 暴 的 观测 温度 和 热能 量 
流量 , 我 们 可 以 计算 出 流出 物 的 初始 半径 mm (从 中 心 引 擎 到 喷 流 喷射 处 的 距离 ) 对 于 ryn > 15 
的 情况 可 以 更 进一步 计算 伽 玛 暴 的 饱和 层 半 径 #， 光 球 层 半 径 mn， 流出 物 滑 行 运动 阶段 洛 伦 
BAF ro MHP 1 式 ， 可 以 判断 光 球 层 和 饱和 层 的 相对 位 置 ， 我 们 发 现 GRB131231A 
中 有 热 成 分 区 间 均 满足 (Fpp/Fwr) 习 /2 > 1, 所 以 我 们 判定 这 个 暴 的 所 有 壳 层 的 光 球 层 位 于 
饱和 层 之 上 (ryn > n), 所 以 我 们 应 用 [26] 中 的 理论 计算 这 个 暴 的 初始 半径 mm 饱和 层 半 径 信 、 
光 球 层 半径 mm、 流出 物 滑行 运动 阶段 洛 伦 兹 因子 ro 

我 们 能 得 到 与 热 成 分 相关 的 参数 值 R, ; R 可 以 解释 为 辐射 区 的 有 效 横 向 尺度 中 ,常数 R, 
意味 着 光 球 的 有 效 辐射 区 域 与 时 间 无 关 。 我 们 可 以 通过 如 下 关系 式 计算 得 到 伽 玛 暴 的 只 值 。 


= FpB 1/2 E (1+2)? Tph 
R= (BEV? = (1.06) HY (4) 


ox 


oa asp 是 斯 特 藩 - 玻 尔 效 曼 常量 ， 单 位 为 erg :.s-1.cm-2 K4; Fpgp 是 观测 到 的 黑体 能 量 流量 , 单 
© fi erg: s-1.cm-?; 7 是 黑体 的 观测 温度 ， 单 位 为 keV，z 为 伽 玛 暴 的 宇宙 学 红 移 ， 必 是 光 
度 距离 ， 单 位 为 cm，r 为 流出 物 滑行 阶段 的 洛 伦 兹 因子 。 

如 图 6 Ca) 所 示 在 脉冲 起 始 阶 段 R 表 现 为 快速 的 上 升 ， 在 22s 后 R 没 有 表现 出 单调 递 
增 的 趋势 ， 变 化 较为 平缓 。 流 出 物 速度 方向 与 视线 的 夹 角 为 6， 在 脉冲 峰值 时 刻 之 后 ，6=0 
区 域 的 辐射 减弱 ， 高 纬 辐射 起 主导 作用 中 ， 此 时 光 球 的 有 效 辐射 区 域 与 时 间 无 关 。 

为 了 分 析 温 度 随 时 间 的 演化 ，Ryde 依据 文章 [28] 的 拐 折 守 律 模型 派生 了 一 种 模拟 温度 
的 拐 折 窜 律 模 型 *。 揭 折 窜 律 的 形式 如 [27] 中 6 式 ， 我 们 使 用 它 去 拟 合 该 暴 探测 到 热 成 分 
区 间 的 温度 与 时 间 , 如 图 6 Cb). 我 们 得 到 对 应 于 6 式 的 参量 拟 合 值 分 别 为 @ = 2.7742. 44、 
to = 24.59 圭 0.95、E= -4.88 士 2. 49、8= 0. 0028 土 1. 08; 由 拟 合 参 量 我 们 得 到 拐 折 容 律 函 
数 的 两 段 上 升 和 下 降 阶 段 寡 律 指 数 分 别 为 a = 7.65 士 3. 48、b = -2. 10 土 3. 48。 
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图 6. a 图 中 展示 RR® 随 脉冲 的 演化 图 ; b 


蓝 色 虚线 为 拟 合 线 。 
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图 中 展示 探测 到 热 成 分 的 光谱 中 温度 随 脉 冲 的 演化 ， 


Fig.6 The figure (a) shows the evolution of R with pulse; The figure (b) shows the evolution of 
temperature with pulse in the spectrum of detected thermal components, and the blue dotted line is 


the fitting line. 


{EAN 385. 6+0. 34,- 


GRB131231A 的 洛 伦 效 因子 呈现 先 上 升 后 随时 间 衰 减 的 演化 ， 起 始 值 为 157 士 0.5， 峰 


它 能 用 拐 折 峰 律 较 好 的 拟 合 ， 如 图 7(a) 上 升 阶 段 的 究 律 指数 u=1. 37 


土 0. 39， 下 降 阶 段 的 暴 律 指数 B.=-1. 38 士 0. 39。S. Iyyani 发 现 具 有 显著 光 球 成 分 的 伽 玛 暴 
中 ， 洛 伦 兹 因子 随时 间 单 调 减 小 r x (F/RD)L4Y14C2P2。 在 脉冲 的 上 升 阶段 ， 总 流量 F 和 R, 
的 增加 速率 几乎 相同 , 洛 伦 兹 因子 接近 固定 值 或 有 时 仅 微 弱 的 随时 间 衰 减 。 我们 所 得 到 的 洛 
伦 效 因子 呈现 为 先 增加 后 衰减 , 不 符合 传统 光 球 模型 预言 的 洛 伦 兹 因子 在 初始 阶段 近似 常数 。 
初始 半径 表示 喷 流 开始 加 速 。 在 GRB131231A 中 no 从 初始 值 (3. 6x10" 士 1. 56x109) 
Y3/2cm 开始 随 光 变 曲线 上 升 ,在 脉冲 峰值 处 (22s ) 达 到 最 大 值 (1. 39x 10°45. 17x10) Y3/2cm， 


但 是 并 没有 很 好 的 追踪 光 变 | 


1 线 ; 在 脉冲 峰值 后 ，m 又 呈现 出 


上 升 的 情况 ， 这 可 能 是 由 于 中 


心 引 擎 的 再 活动 导致 ， 在 脉冲 后 期 m 下 降 (6. 9x10°+1.02x10°) Y3/2zcm， 低 于 初始 值 。m 的 


变化 量 级 分 布 在 10 ~ 10 ， 这 个 演化 范围 在 黑洞 的 视界 半径 CLO" cm; 对 于 一 个 黑洞 质量 5~ 
10Me) 与 Wolf-Rayet 前 身 星 的 核心 大 小 之 间 ， 这 表明 mo 与 喷 流 和 前 身 星之 间 的 相互 作用 有 
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Fig.7 The evolution ofr , ro , 7% and ryn of GRB131231A. 

饱和 层 半径 ns 的 演化 与 初始 半径 在 脉冲 中 的 行为 相似 ， 在 数值 上 相差 接近 三 个 数量 级 。 
初始 值 为 (5. 7x 10°+2. 46x10 Y-5/4cm, x 最 大 值 为 (4. 72x10" 土 1. 76x10”) Y-5/4cm, r 
演化 过 程 中 的 最 小 值 为 (1. 2x10°+1. 83x10°) Y-5/4cm。 

光 球 层 是 可 以 观察 到 光子 的 最 深 区 域 。 在 脉冲 起 始 至 达到 峰值 阶段 (22s )， 光 球 层 半 
径 rpn 呈 现 跟随 脉冲 上 升 ( 追 踪 )， 从 初始 值 (2.06x10" 土 8.68x10") Y-14cm 在 近似 脉冲 峰 
值 处 达到 最 大 值 (1.1x10” 圭 5.89x10") Y-14cm; 脉冲 峰值 后 , 7 快速 下 降 ， 变 化 范围 在 一 
个 量 级 内 ， 之 后 仅 微弱 变化 。GRB131231A 的 光 球 层 ， 最 小 值 为 (1.41x10" +6.87x10") 
Y-14cm。 脉冲 下 降 阶 段 rn 变化 与 GRB110721A °”! GRB140329B 和 GRB160625B "一 致 。 
rph K SKS, 7 为 伽 玛 暴 总 的 运动 学 光度 ，M 是 重子 载 量 ,满足 M = L/T; rpn 的 变化 与 
洛 伦 兹 因子 [和 M 有 关 。 脉 冲 的 上 升 阶段 ，mpn 随 脉冲 上 升 ，[ 也 随时 间 增 加 ， 这 要 求 Lo 或 M 要 
以 更 快 的 速率 增加 ， 而 重子 载 量 的 增加 ， 导 致 与 它 相关 的 电子 增加 了 不 透明 度 ， 从 而 增加 了 
光 球 层 半径 。 在 脉冲 衰减 阶段 ，[ 和 Lo 都 减少 ， 而 mn 保持 常数 或 仅 微 弱 地 变化 ， 为 了 出 现 这 
种 状况 ， 就 要 满足 Lo x T3 (或 M xT?) PI, 


os 


4. 结 论 与 讨论 


GRB 131231A 是 一 个 比较 明亮 且 较 为 接近 单 脉冲 的 长 暴 。 我 们 应 用 了 贝 叶 斯 方法 ， 
利用 经 验光 谱 模型 和 物理 模型 ， 即 CPL 模型 、Band 模型 和 普 朗 克 函 数 模型 (BB )， 执 行 详 
细 的 时 间 分 辨 谱 分 析 , 并 在 每 个 时 间 分 辨 谱 区 间 上 进行 模型 比较 , 来 确定 适合 某 一 区 间 的 最 
优 模型 ， 我 们 总 共 获 得 了 24 个 光谱 并 研究 了 它们 的 性 质 。 这 个 伽 玛 暴 的 所 和 wa 值 都 表现 为 
较 好 的 追踪 行为 。 本 的 流量 追踪 行为 可 能 与 两 种 同步 辐射 模型 有 联系 [571 。 图 2(b) 中 显示 
a 的 值 均 分 布 于 同步 辐射 死 线 之 下 ， 整 体 分 布 非常 符合 同步 辐射 模型 ， 李 亮 等 人 指出 在 上 升 
阶段 ，a 变 得 更 硬 ， 这 表明 电子 正在 从 快速 冷却 区 向 慢 速 冷却 区 演化 四。 在 衰减 的 阶段 ，a 
在 观察 上 变 得 更 软 ， 宽 脉冲 的 衰减 阶段 可 能 由 “曲率 效应 ”控制 CIB303]。 
通过 最 优 模型 与 BB 模型 的 共同 拟 合 , 我 们 探测 到 16 个 有 明显 热 成 分 的 时 间 分 辨 谱 区 
间 ， 该 暴 不 是 热 主 导 的 伽 玛 暴 。 李 亮 系统 研究 了 8 个 含有 热 成 分 伽 玛 暴 的 能 谱 特 性 站， 对 于 
非 热 成 分 主导 的 伽 玛 暴 ， 他 发 现 加 入 了 热 成 分 后 , a 值 变 软 ，Ee 变 硬 ， 能 量 流量 F, 在 不 同 的 
模型 下 没有 明显 的 不 同 。 对 于 GRB 131231A， 在 光 变 曲线 峰值 之 前 ， 加 入 BB 成 分 后 Ep 值 
变 硬 ; 峰值 之 后 ， 加 入 BB 成 分 后 Ep 值 变 软 。 在 光 变 曲 线 峰值 之 前 ， 加 入 BB 成 分 后 a 值 变 
软 ， 峰 值 之 后 ， 加 入 BB 成 分 后 变 硬 。GRB131231A 具有 特殊 的 光谱 演化 。 

GRB 喷 流 的 中 央 引 警 很 可 能 是 一 个 混合 系统 ， 既 有 一 个 热 的 火球 成 分 ， 又 有 一 个 冷 的 
坡 印 廷 流 成 分 史 。 一 些 明亮 的 暴 〈 例 如 GRB 090902B ) 在 脉冲 的 开始 处 或 峰值 附近 表现 出 
类 似 热 的 时 间 分 辨 谱 , 在 随后 的 脉冲 中 出 现 非 热 的 时 间 分 辨 谱 , 这 可 能 表明 喷 流 成 分 的 逐渐 
AFALE, 高 稚 等 人 在 解释 伽 玛 暴 的 喷 流 成 分 的 起 源 方面 做 了 更 全 面 的 描述 ， 即 混合 喷 流 模 
型 60。 李 亮 使 用 此 方法 分 析 了 8 个 Fermi GRB! ， 发 现 大 多 数 爆发 可 以 用 混合 喷 流 情况 很 
好 地 解释 。GRB 131231A 的 低能 光子 指数 a 的 值 不 超过 所 谓 的 同步 辐射 死 线 ， 但 加 入 热 
成 分 后 ， 拟 合 统计 量 明显 改善 ， 其 中 ABIC <-8 的 光谱 主要 集中 脉冲 开始 阶段 和 峰值 阶段 
(10~26s )， 随 后 热 成 分 逐渐 减少 ， 说 明 伯 玛 暴 可 能 从 物质 主导 火球 过 渡 到 磁 主导 的 坡 印 廷 
流量 (Poynting-Flux);， 在 脉冲 后 期 又 出 现 热 成 分 ， 这 可 能 是 中 心 引擎 再 活动 导致 。 
虽然 温度 也 能 用 拐 折 寡 律 拟 合 , 但 是 与 Ryde 拟 合 伽 玛 暴 中 热 脉冲 所 得 到 的 上 升 阶段 和 
下 降 阶 段 的 指数 值 均 存 在 较 大 差异 C。 

我 们 计算 了 与 热 成 分 相关 参数 r、ro、7s、rpn。 我 们 得 到 的 洛 伦 兹 因子 呈现 为 先 增加 后 


A dh 


EHA 


FR 


Tij 
J 


衰减 ， 不 符合 传统 光 球 模型 预言 的 洛 伦 兹 因子 在 初始 阶段 近似 常数 请 P2， 之 后 以 指数 衰减 
的 情况 ， 这 可 能 是 因为 内 激 波 模型 导致 50， 另 一 种 解释 是 孤立 磁 星 诞生 的 情况 ， 此 时 洛 伦 
兹 因子 随时 间 增 加 0 但 是 GRB131231A 的 光度 远大 于 磁 星 模型 的 光度 , 所 以 磁 星 模型 的 设 
想 不 太 可 能 成 立 。 我 们 得 到 mo、 和 wy 的 参数 值 的 变化 范围 与 具有 显著 光 球 成 分 的 伽 玛 暴 
GRB 081224A、GRB 090719A、GRB 100707A 等 22P9 变 化 范围 吻合 。 其 中 ，7m 的 演化 范围 


在 黑洞 的 视界 半径 (10° cm 对 于 一 个 黑洞 质量 5~ 10Me) 与 Wolf-Rayet 前 身 星 的 核心 大 小 


之 间 , 这 表明 mo 与 喷 流 和 前 身 星之 间 的 相互 作用 有 关中 9。 ns 的 演化 行为 与 近似 , 在 数值 上 相 
差 三 个 数量 级 之 内 。 在 脉冲 起 始 至 达到 峰值 阶段 (22s )， 光 球 层 半径 my 呈现 跟随 脉冲 上 升 
(追踪 ); 脉冲 峰值 后 ，7pn 快 速 下 降 ， 变 化 范围 在 一 个 量 级 内 ， 之 后 仅 微弱 变化 。 脉 冲 的 上 
升 阶段 ，7ppn 随 脉冲 上 升 ，[ 也 随时 间 增 加 ， 这 要 求 Lo 或 M 要 以 更 快 的 速率 增加 ， 而 重子 载 量 
的 增加 ， 导 致 与 它 相关 的 电子 增加 了 不 透明 度 ， 从 而 增加 了 光 球 层 半 径 。 在 脉冲 衰减 阶段 ， 
rT 和 Lo 都 减 小 ， 而 rjn 保 持 常数 或 仅 微弱 地 变化 ， 这 种 状况 就 要 求 Lo x T3 (或 M x T2 ) P, 

GRB 131231A 的 能 谱 特 征 行为 表明 这 个 暴 有 一 些 特殊 性 ， 需 要 寻求 合适 的 
理论 进行 解释 ， 值 得 进一步 深入 研究 。 
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Study on energy spectrum characteristics of Ep and a 


"double tracking" gamma-ray bursts 


Wang Chen, Peng Zhao-Yang, Chen Jia-Ming 
(College of Physics and Electronic Information, Yunnan Normal University, Kunming 650500, China, 
Email: pengzhaoyang412@ 163.com) 

Abstract: GRB 131231A is a special gamma-ray burst (GRB) with Ep and a having "dual 
tracking" behavior. To further analyze the energy spectrum characteristics of the burst, we used the 
Bayesian method to perform a detailed time-resolved spectral analysis and found that 16 of the 24 
resolved spectra have significant thermal components, which are mainly distributed in the early 
and peak stages of the pulse. Near the peak, the transition from matter-dominated to 
magnetic-dominated GRB jets is supported. After adding thermal components, parameters Ep and 
a become hard, while the peak flux F, is almost unchanged. We calculated the characteristic 
parameters of the relativistic flow using the composition of the photosphere, and found that the 
initial radius ro of the flow characteristic parameter, the radius of the saturation layer Ts, the radius 
of the photosphere layer rpn, and the Lorentz factor r during the coast phase of the flow all 
occurred significantly. The evolution of temperature over time can also be fitted by the inflection 
power-law model. The power-law index of its rising and falling stages are a = 7.65+43.48 and b = 
-2.10+3.48, respectively. But b has a large deviation from the index value close to -2/3 obtained 
by previous study of other thermal pulses. 


Key words: Fermi Gamma-ray burst; Prompt emission; Spectral fitting; Data analysis; optimal model 


表 4 伽 玛 暴 131 231A 最 优 模 型 拟 合 的 时 间 分 辨 谱 
Table.4 The best fitting results of models to fit the time-resolved spectrum of GRB 131231A 


Toar Toni Models S/N a Ep(keV) kT(keV) Fp(erg s'cm’) Fpp(erg scm) 

12.0: 14.17 CPL’ 29 -0.95%893  162.28+343 7.46E — 07+375E-98 
CPL7+BB —1.02+9:11 196.44372 21.32+238 = -7.62E — 07243408 2.03E — 0874-42808 

4.17:15.32 CPL? 31 一 1.05+0.07 284.02+52.0 1.36E — 06+726E-08 
CPL’+BB —1.09+0.08 329.6+264 22.721631 = 1.37E — 0677582708 3.07E — 0875-678-08 

15.32:18.03 CPL? 85 —1.0+0.03 346.34+25.0 3.29E 一 06+612E-08 
CPL+BB —1.07+0.03 492.7+668 24.20+215 3.38E— 06+5.68E-08 2.10E — 07+489E-08 

18.03:20.08 Band 108  —0.77+0.98 207.3414 4.93E — 0613328708 


Band + BB —1.0479-96 388.3+30245 25.344187 5.08E — 06+9.82E-08 4.62E — 07+1.95E-07 
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3192 
3184 
2331 
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20.08:21.35 
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40.97: 
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31.01 


32.61 


33.45 


35.41 


37.26 


38.65 


39.61 
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42.01 


45.0 
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Band + BB 
CPL? 
CPL2+BB 
CPL? 
CPL2+BB 
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Band 
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Band 
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Band 
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+ BB 
Band 


Band 
+ BB 
Band 


Band + BB 
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85 


94 


159 


83 
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165 
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144 


91 
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94 


64 


64 


44 


58 


—0.83+8.84 
-0.814943 
一 0.93+8.83 
—0.97 7598 
一 0.91+8.83 
—0.93+8.83 
—0.87+8.82 
—0.9175-93 
—0.86+8.84 
—0.63+933 
—0.81+8.84 
—0.69+8.88 
一 0.77+8.84 
—0.61+8.84 
—0.83+8.08 
—0.72+8.87 


—0.84+8.04 


—0.84+8.05 


—0.89+8.84 


—0.8525:88 


—1.08+8.98 


一 1.01+6.33 
—1.08+8.98 
—1.01+8.03 
一 1.29+6.04 
—1.28+6.06 
一 1.36+6.04 
—1.13+933 


一 1.34+8.8 


—1.274333 
—1.3723-08 
-1.274 
—1.44+8.08 
—1.44+8.86 
—1.50+8.08 
—1.325338 
—1.56+8.08 
—1.534333 
—1.53+8.08 


1 5+0.07 
1.57508 


233.60+17:9 
234.41059 
360.864369 
481.0+64.5 
389.18+269 
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375.444228, 
220.877189 
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1.48E 一 05+ 一 07 
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9.09E — 06+121E-07 
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1.08E — 05+131E-07 
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7.99E — 06+127E-07 
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2.35E — 06+321E-08 
2.38E — 06+447E-08 
1.88E — 06+401E-08 
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2.55E — 064395708, 
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3.91E — 07+2.07E-07 


2.37E — 08+1.04E-07 


2.87E — 07+115E-07 


1.81E — 08+2.36E-07 


9.81E — 09+8.57E-08 


2934 
2912 
1811 
1777 
1754 
1736 
3014 
3003 
1188 
1171 
2480 
2462 
3315 
3266 
2062 
2044 
2281 


2275 


3529 
3521 


3142 
3141 


3193 
3194 
2025 
2019 
3373 
3351 
3091 
3095 


2767 
2755 
2116 
2108 
2525 
2512 
2056 
2046 
3647 


3638 
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